Transfer reactions and the dispersive optical-model 
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The dispersive optical-model is applied to transfer reactions. A systematic study of (d, p) reac- 
tions on closed-shell nuclei using the finite-range adiabatic reaction model is performed at several 
beam energies and results are compared to data as well as to predictions using a standard global 
optical-potential. Overall, we find that the dispersive optical-model is able to describe the angular 
distributions as well as or better than the global parameterization. In addition, it also constrains the 
overlap function. Spectroscopic factors extracted using the dispersive optical-model are generally 
lower than those using standard parameters, exhibit a reduced dependence on beam energy, and are 
more in line with results obtained from (e, e'p) measurements. 



I. INTRODUCTION 

Rare isotopes have posed new challenges to our un- 
derstanding of nuclei as protons and neutrons organized 
into well-defined orbitals. Nucleon-nucleon (NN) cor- 
relations manifest themselves in different forms as one 
moves away from the valley of stability. Nuclear reac- 
tions form one of the most important probes for the lim- 
its of stability and the role of the underlying fundamental 
interactions. Transfer reactions offer a plethora of op- 
portunities for studying shell structure because one can 
access ground as well as excited states (including res- 
onances), and by choosing appropriately the kinematic 
conditions, one can explore peripheral properties as well 
as the nuclear surface. For this reason, they consti- 
tute an important part of the science program of many 
rare isotopes facilities worldwide. Examples of recent 
studies include 132 Sn(d,p) 133 Sn Q, 15 C(d,p) 16 C 0, and 

34,36,46 Ar(p,d) @. 

The increased interest in using transfer reactions as a 
tool to study these exotic nuclei has called for reaction- 
theory developments [H-Il]. Whereas in the sixties the 
standard approach to one-nucleon transfer was to rely 
on the first-order distorted-wave Born approximation 
(DWBA), it is now well understood that in (d,p) re- 
actions, breakup is important and first-order perturba- 
tion is insufficient to describe the process. A practical 
method for including deuteron breakup was introduced 
by Johnson and Soper 0, within a zero-range approx- 
imation, usually referred to as the adiabatic wave ap- 
proximation ( ADWA) . In addition to including deuteron 
breakup to all orders, the ADWA depends exclusively 
on nucleon optical-potentials rather than the more am- 
biguous deuteron optical-potential as in DWBA. ADWA 
was later extended by Tandy and Johnson [ljj to in- 
clude finite-range effects (FR-ADWA). A recent system- 
atic study of [d,p) reactions within FR-ADWA [12| has 
shown the importance of finite-range effects in consid- 
ering deuteron breakup in (d,p) reactions. The method 
explored in Refs. [ill EH is based on the truncation of a 
Weinberg expansion. In Ref. [l3j a systematic compar- 



ison between (d,p) angular distributions for FR-ADWA 
and those from the exact full Faddeev solution is per- 
formed. The results from FR-ADWA [U| are within 10% 
of the full solution at forward angle, for a wide range of 
beam-energies 13J. While the exact full solution is com- 
putationally intensive and expensive, FR-ADWA calcu- 
lations are clearly of practical use. 

The work in [l]| determines the level of accuracy of 
FR-ADWA concerning the solution of the three-body 
n + p + A dynamics in the reaction. For the (d,p) re- 
actions on 48 Ca at 19 MeV and 56 MeV effects are 6% 
and 3% respectively [l3|. With this level of accuracy it 
is now essential to have good control over the ingredients 
to the problem, namely the effective interactions. All 
applications of FR-ADWA so far have relied on global 
optical potentials for the nucleon-target interactions and 
"standard single-particle parameters" for the effective in- 
teraction that determines the final neutron bound state. 
So far, these two types of input have been disconnected. 
Future extraction of nucleon properties in rare isotopes 
will require presently unavailable knowledge of nucleon- 
target interactions in the continuum as well as bound- 
state information. 

An attractive feature of the recently implemented dis- 
persive optical-model (DOM) [14Ul6l ] is that it intrin- 
sically connects the scattering and bound states. The 
method, first introduced by Mahaux and Sartor fl7| , and 
also known as the dynamic polarization potential (DPP), 
relies on the dispersion relation between the imaginary 
and real parts of the nucleon self-energy. The work of 
Ref. [HI concentrated on the interaction between pro- 
tons and Ca isotopes and found larger correlations near 
the Fermi surface with increasing nucleon asymmetry. 
Neutrons were included in the fits reported in Ref. [la ]. 
That analysis raised questions concerning the standard 
parameterization of the imaginary part of the nucleon 
potentials as a function of asymmetry {N — Z)/A, par- 
ticularly for neutrons. More recently, a non-local exten- 
sion of the DOM has been introduced that allows for a 
broader application to calculate observables below the 
Fermi energy [l8j . Insight from ah initio calculations of 
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the nucleon self-energy are expected to provide further 
guidance on the details of non- locality jl9j |. 

In the present work we will use only the local version 
of the DOM. Its most recent implementation [l6| yields 
an accurate representation of a wealth of data in several 
locations of the chart of nuclides, including nuclei with 
Z = 20,28 and N = 28 (FIT1), nuclei with N = 50, 
with Z = 50 (FIT2), and finally Z = 82 (FIT3). When 
several isotopes or isotones are included in the fit, the 
observed nucleon asymmetry dependence of the poten- 
tials allows for an extrapolation to nuclei that are more 
exotic [l6j . The corresponding predictions can thus be 
probed with future experiments that may allow further 
refinements. This feature of data-driven extrapolations 
to the respective drip lines and its capacity to link both 
reaction and structure information make the DOM an 
excellent framework to study rare isotopes. 

In DOM fits [l4l - [l6j . nucleon elastic, total and reac- 
tion cross section data, as well as bound-state properties 
inferred from (e, e'p) measurements, are all included si- 
multaneously. The aim of the present work is to explore 
the application of the DOM to transfer reactions. We in- 
clude a number of representative cases: 40 Ca((i,p) 41 Ca, 
48 Ca(d,p) 49 Ca, 132 Sn(d,p) 133 Sn, and 208 Pb(d,p) 209 Pb. 
These reactions are studied at various beam energies for 
which data are available. We note that apart from the 
n — p interaction that binds the deuteron, all effective 
interactions to the problem in the ADWA framework 
are provided by the DOM. We compare the results to 
those obtained using a global optical-potential (CH89 
from Rcf. (2Cj) and standard geometry for the bound 
state. 

In Sec. [II] we summarize the important theoretical con- 
cepts pertaining to the ADWA and DOM. Details of the 
ingredients of the reaction description are presented in 
Sec. IIII Al and the results of the transfer calculations are 
presented in Sec. IIII Bl A more detailed discussion on 
spectroscopic factors is presented in Sec. IIII CI Finally, 
we summarize and draw our conclusions in Sec. [IVJ 



II. THEORY 

A. Finite-range adiabatic-wave approximation 

The adiabatic theory of Rcfs. 0,[ll| for A(d,p)B starts 
from a three-body model of n+p + A. The deuteron scat- 
tering wavefunction in the incident channel is obtained 
by solving the differential equation: 

[E + ie - T r - T R - U nA - U pA - V np ]& + \r, R) 

= ie<j) d (r) cxp(iK d ■ R), (1) 

with r = r p — r n (R = (r n + r p )/2) being the rela- 
tive coordinate (center-of-mass coordinate) of the n — p 
system. The neutron and proton coordinates, which are 
taken at the center of mass of the target A, are given 
by r n and r p , respectively. We take U nA (r n ), U pA (r p ), 



and V np (r) to be the neutron-target, proton-target, and 
neutron-proton interactions. In this work, we use the 
Reid potential [2l| for V np and DOM potentials for U pA 
and U nA (see Sect. IIIBj) . 

In this three-body approach, the solution of Eq. ([!} is 
inserted into the exact transfer matrix element: 

T = (<t>nAX ( pB } I V np + A rem |tf <+>) , (2) 

where <p nA is the bound state of the neutron-target sys- 
tem, and x p b 1S a P ro ton scattering distorted wave in 
the outgoing channel. The remnant operator is A rem = 
U p a ~ U p b- Contributions from this term are often small 
except for the lighter systems. 

Johnson and Tandy [ll| noted from Eq. @ that the 
exact three-body wavefunction \]>( + ) (r, R) is only needed 
within the range of the V np interaction, whenever the 
remnant contributions are negligible. For this reason, 
the Weinberg states, which form a complete square in- 
tegrable basis within the range of V np , offer an excellent 
representation for the three-body scattering wavefunc- 
tion: 

oo 

¥+\r,R)=Y J Hr)^{R)- (3) 

i=l 

Here 4>i(r) denote Weinberg states, normalized by 
(<^i|Kip|0j) — — &ij ancl satisfying: 

(T r + a t V np )Mr) = -td<t>i{r), (4) 

where e d is the deuteron binding energy and a, are the 
eigenvalues. Inserting Eq. (|3|) into Eq. ([T]) one arrives at 
a non-trivial coupled-channels scattering equation. How- 
ever, if only the first term in the Weinberg expansion 
of Eq.® is necessary to describe transfer, then Eq.([T]) 
becomes a single-channel optical-model-type equation: 

{E + e d + ie-T R -Uj T (R))x{ T (R) = ieN d exp(iK d -R), 

(5) 

where N d is the normalization coefficient defined as N d = 
~{4>i\V np \4>d)- The potential Ujt in Eq.© is calculated 
by a folding procedure: 

Ujt(R) = -(Mr)\V np (U nA + U pA )\Mr)) (6) 

where, 4>i{r) is the lowest eigenfunction of Eq. Q, pro- 
portional to the deuteron bound-state wavefunction 0^. 
The expansion of Eq.Q truncated to the first term, is 
now inserted into the transfer matrix element : 

T = (<t>nAX p ~J\V np + A rem \Mx( T (R)/N d }) . (7) 

The validity of the truncation has been tested through 
comparisons with Faddeev solutions and the results sug- 
gest that for most cases of interest, the finite-range 
ADWA works very well [H] . 

One important advantage of ADWA over DWBA, is 
that ADWA depends on nucleon optical-potentials rather 
than the deuteron optical-potential, which is far more 
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ambiguous. ADWA also depends on the nuclcon bind- 
ing potential, as does the DWBA. Since these nucleon 
optical and binding potentials are connected by a dis- 
persion relation in the DOM approach, it is interesting 
to test how the DOM potentials perform in describing 
transfer angular distributions for a wide range of ener- 
gies and targets. In the next subsection we give a brief 
summary of the recent developments with DOM nucleon 
optical-potentials. 



and a similarly scaled dispersive correction. For lack of 
detailed theoretical input, the imaginary part of the self- 
energy is also approximated by a local potential in keep- 
ing with a long-standing tradition of empirical optical 
potentials. We note that the dispersive correction is cor- 
respondingly local. The nucleon self-energy, as typically 
represented in the DOM, can then be written as 

U(r,E) = V H f(t,E) + AV(r,E) +iW(r,E). (12) 



B. Dispersive optical-model (DOM) 

The dispersive optical-model [l7| develops a practical 
representation of the irreducible nucleon self-energy E* 
from the Green's-function approach to the many-body 
problem [22[. The real part of the nuclcon self-energy 
or optical-model potential can be decomposed into an 
energy-independent non-local part and a dynamic con- 
tribution (energy-dependent), that can also be non-local, 



Re S (r, v'\ E) = Re £ (r, r';e F ) + AV(r, r'\ E), (8) 

where Ef is the Fermi energy. The second term, referred 
to as the dispersive correction, can be obtained from the 
imaginary part through the subtracted dispersion rela- 
tion 



AV(r,r';£) = 

+ -V J Im S (r, r'; E 1 ) 



(9) 



E' — E E'-e F 



dE', 



where V stands for the principal value and we note the 
convention to employ the same sign for the imaginary 
part of the self-energy above and below the Fermi en- 
ergy 17 1 . By definition in Eq. ([8]). the dispersive correc- 
tion is zero at the Fermi energy. The dispersive correction 
varies rapidly around Ep and causes the valence single- 
particle levels to be focused towards the Fermi energy. 

Following a long tradition [23j], the non-local energy- 
independent term ReS(r, r'; Ef) is approximated by a 
local energy-dependent form, generally designated as the 
Hartree-Fock (HF) potential VW(r, E 1 ). While it does 
describe mean-field properties, it is not strictly a HF 
contribution, since it involves the complete self-energy 
at one particular energy. The energy derivative of Vhf 
encodes a measure of the non-locality, which is related to 
the momentum-dependent effective mass 



m (r, E) 



= 1 - 



dV HF (r,E) 
dE 



(10) 



where m is the nucleon mass. 

The local approximation of the HF potential necessi- 
tates a scaled imaginary potential given by 



w= m(r,E) ImE 



m 



(11) 



The Fermi energy is defined as 



e F 



2 

M A+1 - (M A + m) 
M A -(M A -i + m), 



(13) 

(14) 
(15) 



where and ej. represent the binding energy for adding 
or removing a nucleon, or alternatively, the single-particle 
energies of the valence particle and hole states. 

Since the current implementation of the DOM includes 
scattering data up to 200 MeV, a lowest-order relativistic 
correction is employed in solving the radial wave equa- 
tion m 



d?_ 
dp 2 



i - 



U(p,E) 



£(£ + 1) 



M 



u(p) = 
(16) 



with p = kr, where k = j^y/T (T + 2m), T is the 
laboratory kinetic energy, E to t is the total energy in the 
ccnter-of-mass frame, and M is the target mass. The 
scaled potential is 



U = 7^,7 



2 (Etot ~ M) 
E tot -M + m 



(17) 



We denote by (pnij ( r ) bound-state solutions to the ra- 
dial wave equation, where n refers to the corresponding 
state in the yl±l-system. The normalized wave functions 
corrected for non-locality are then given by 



<Pnij ( r ) 



m(r,£ nlj ) _ 
— <PrUj (r) ■ 



(18) 



where e n ij is the discrete energy solution to Eq. (|T6"|) . It 
should be emphasized that these corrected wave functions 
correspond to overlap functions for adding or removing 
particle from the target nucleus [l8j], e.g. 



\/Snej<Pnej(r) 



l*0> 



(19) 



where ip n £j is normalized to one and S n £j represents the 
spectroscopic factor (norm of the overlap function) . The 
solutions Tp therefore still need to be normalized by the 
spectroscopic factor which is obtained in the DOM from: 



dr. 



(20) 
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where the energy-dependent effective mass determines 
the reduction of strength from the mean-field value and 
is given by 



m (r, E) 



1 



dAV(r, E) 



m (r, E) 



dE 



(21) 



The result for the spectroscopic factor in Eq. (|20[) was 
derived in Ref. [llj and was shown to be an excellent 
quantitative approximation to the corresponding solution 
of the Dyson equation that incorporates a non-local HF 



in Ref. |18|. We will also include results for the root- 



mean-square (rms) radius given by 



T>rms 



Vllj 0) r 2 dr 



(22) 



in the DOM analysis. 

In the present work, we will extract the normalization 
of the overlap functions from the analysis of the transfer 
reactions by scaling the theoretical cross section to match 
the data. We will also compare these results with the 
DOM spectroscopic factors generated by Eq. (|20[) . For 
comparison purposes, we will also consider the single- 
particle wave function produced with a Woods Saxon 
potential with standard parameters (radius ro = 1.25 
fm and diffuscncss a = 0.65 fm) and the depth adjusted 
to reproduce the experimental binding energies. These 
wave functions will be referred to as <pws- An important 
additional experimental constraint on the overlap func- 
tion at large distances (when available) is provided by 
an asymptotic normalization constant (ANC). The sin- 
gle particle ANC b n £j is defined through the asymptotic 
behavior of the single particle state ip 



(23) 



where K = ^2/j,s n ij , fi is the corresponding reduced mass 
and tp n ej(f) is mormalized to one. We note that the 
spherical Hankcl function in Eq. (|23p is only appropriate 
for neutrons. 



III. RESULTS 

A. Details of the calculations 

Separate DOM fits were produced for different parts 
of the chart of nuclides in Ref. [HI]. For the Ca isotopes, 
the DOM fit included nucleon scattering and bound-state 
data for nuclei with Z = 20,28 and N = 28 (FIT1). 
Excluded were proton clastic scattering data on 40 Ca 
below 18 McV, since for those cases the optical model 
does poorly and there are fluctuations in the reaction 
cross section [2]|. For the Sn isotopes, the DOM fit in- 
cludes proton clastic-scattering data on 112_124 Sn and 
neutron elastic-scattering data on 116 Sn, 118 Sn, 120 Sn, 
and 124 Sn (FIT2) Finally, for Z = 82, 208 Pb data 



TABLE I. Properties of overlap functions with a compari- 
son of the DOM overlap function, corrected for non-locality 
(p with the Woods-Saxon single-particle wavefunction tpws- 
These include the counting number n, the angular momenta 
(Ij) of the valence orbital, the separation energy S„, the root 
mean square radius of the valence orbital R rms , and the mod- 
ulus of the single-particle ANC b n ij . 



Nucleus Overlap 


n£j 


S n [MeV] 


R rms [fm] 


[fm- 1/2 ] 


41 Ca 


Vws 
V 


O/7/2 


8.362 


3.985 
3.965 


2.285 
2.261 


49 Ca 


Vws 
V 


lf>3/2 


5.146 


4.606 
4.820 


5.818 
6.098 


133 Sn 


Vws 
V 


1/7/2 


2.469 


6.080 
6.513 


0.844 
0.831 


209p b 


Vws 
V 


lffg/2 


3.936 


6.498 
6.746 


1.650 
1.827 



are included (FIT3). The fitting procedure follows pre- 
vious works [IJ, [l5| but the depth of the HF potential 
was adjusted to reproduce the valence particle levels for 
neutrons. Since transfer cross sections are strongly de- 
pendent on the binding energy of the final bound state, 
particular effort was made in all these fits, to reproduce 
the experimental binding energies. The DOM parame- 
terizations are subsequently used to calculate A{d 1 p)B 
within FR-ADWA. The DOM thus provides all optical 
potentials U n A, U p a, U p b as well as the neutron bound- 
state interaction V n A- The properties of the neutron 
states considered in this work are summarized in Table U 
including the rms radius and the single-particle asymp- 
totic normalization coefficient [1] [see Eq. (f23))]. We use 
TWOFNR [2BI to calculate the finite-range deuteron adia- 
batic potential and FRESCO [27| to calculate the transfer 
cross sections. 

The DOM fit for Ca isotopes generated a radius param- 
eter of 1.18 fm while the corresponding standard Woods- 
Saxon potential used for comparison was fixed at 1.25 
fm. The non-locality correction almost completely can- 
cels this difference making the DOM overlap function 
(p for 41 Ca essentially identical to <pws- When a node 
is present, the non- locality correction is even more pro- 
nounced making the rms radius of the IP3/2 DOM wave 
function larger than its Woods-Saxon counterpart. We il- 
lustrate this in FigO]for the 49 Ca ground state, where the 
square of the single-neutron overlap function obtained 
with the local DOM, corrected for non-locality, is com- 
pared to the Woods-Saxon counterpart with standard ge- 
ometry. The DOM result (FIT2) for the radius parame- 
ter for Sn nuclei is 1.24 fm and therefore a substantially 
larger rms radius is obtained for the DOM overlap func- 
tion compared to the standard Woods-Saxon one when 
the non-locality correction is applied for the I/7/2 orbit 
(see Table H|. 

As is usual in optical-model potentials, the imaginary 
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FIG. 1. (Color online) Comparison of the square of the single- 
neutron overlap functions for 49 Ca obtained with a Woods- 
Saxon potential ipws (solid) and with the local DOM cor- 
rected for non-locality ip (dashed). 



potential in the DOM is separated into surface and vol- 
ume terms. The former is stronger at the lower energies 
considered in this work and the latter becomes dominant 
only at much higher energies. The volume term was as- 
sumed to have a small dependence on ±(N—Z)/A, where 
the positive and negative signs refer to protons and neu- 
trons. The magnitude of this term was determined from 
the fit to the lead data (FIT3 of Ref. 0) and then was 
kept fixed when fitting the Sn and Ca data [l6|. It is for 
the study of n- 132 Sn (derived from FIT2 in Ref. [H|]) 
that the asymmetry dependence becomes critical, as one 
is extrapolating well outside the region of fitted data. As 
previously seen for the Ca isotopes [15|, and for FIT1 
in Ref. [ly], the Sn DOM potentials contain an insignif- 
icant asymmetry dependence of the neutron imaginary 
surface term, while that for protons exhibited a linear 
dependence on (N — Z)/A. 

A comparison between the DOM and CH89 optical 
potential reveals some similarities and some differences. 
Since the surface absorption for neutrons exhibits a dif- 
ferent (weaker) asymmetry dependence than assumed in 
Ref. [201, it is not surprising that the 132 Sn potentials 
show the largest differences, as this corresponds to a 
substantial extrapolation from potentials constrained by 
data (see Fig(2]). In particular, the neutron (proton) sur- 
face absorption in the DOM is substantially less (more) 
than in the CH89 parameterization). The real parts of 
the CH89 potentials have the same radius parameter for 
both protons and neutrons by decree, whereas in the 
DOM, due to the different surface absorption, the dis- 
persive correction makes the real proton potentials ex- 
tend farther than those for neutrons. We also note that 
the CH89 potentials are not dispersive, making a more 
detailed comparison less productive. 



CD 

> 




r[fm] 



FIG. 2. (Color online) Optical potentials for n- 132 Sn and p- 
132 Sn (thin) at 4.7 MeV: comparing DOM with CH89. The 
real component is shown in panel (a) and the imaginary part 
in panel (b). 



B. Transfer cross sections 



We consider 9 reactions for which data exists: 
40 Ca(d,p) 41 Ca at E d = 20 and 56 MeV, 48 Ca(ii,p) 49 Ca 
at E d = 2,13,19.3 and 56 MeV, 132 Sn(d,p) 133 Sn at 
E d = 9.46 MeV, and 208 Pb(<i,p) 209 Pb at E d = 8 and 
20 MeV. 

FR-ADWA calculations for all these cases were per- 
formed for three interaction models: optical potentials 
from CH89 (20| and the neutron overlap function as pws 
(CH89+WS); optical potentials from DOM but the neu- 
tron overlap function as ipws (DOM+WS); and finally 
both the optical potentials and the neutron overlap func- 
tion from this local approximation of DOM corrected for 
non-locality (p (DOM). In Figs. [3l [4j and [5] the angular 
distributions for E d = 2, 19.3 and 56 MeV, normalized to 
data, are plotted for the 48 Ca(<i, p) 49 Ca reaction. 

The normalization is performed at the peak of the data, 
i.e. at backward angles for reactions below the Coulomb 
barrier and at the first peak when the bombarding en- 
ergy is above it. For these cases the relevant angles were 
9 w 155°, 9 w 10° and 6 « 5° for E d = 2, 19.3 and 56 
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FIG. 3. (Color online) Angular distributions are shown for the 
reaction 48 Ca(d,p) 49 Ca at Ed = 2 MeV. Theory predictions 
have been normalized to the data at backward angles. 



FIG. 6. (Color online) Angular distributions for the 
1,!2 Sn(d,p) 133 Sn reaction at a deuteron energy of Ed = 9.46 
MeV are shown normalized at the peak of the experimental 
cross section. 



40 




6 [deg] 

FIG. 4. (Color online) Angular distributions are displayed 
for the reaction 48 Ca(d,p) 49 Ca at E d = 19.3 MeV. Theory 
predictions have been normalized to the data at the peak. 
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e [deg] 

FIG. 7. (Color online) Angular distributions for the 
208 Pb(d, p) 209 Pb reaction are shown at a deuteron energy of 
Ed = 20 MeV and normalized at the peak of the experimental 
data. 




e [deg] 

FIG. 5. (Color online) Angular distributions for the reaction 
48 Ca(d,p) 49 Ca at E d = 56 MeV are displayed. Theory pre- 
dictions have been normalized to the data at forward angles. 



MeV, respectively. Note that at 2 MeV, compound con- 
tributions to the transfer cross section can be important. 
However for this close analysis of the experimen- 

tal angular distribution suggests such a contribution to 
be negligible. 

The first thing to note is that DOM is able to describe 
the cross sections well. However, there is no significant 
difference between the angular dependence predicted by 
DOM, DOM+WS and CH89+WS, even though there are 
differences in the shapes of the overlap functions as shown 
in Fig. [TJ As a similar set of data was used to constrain 
the CH89 optical potential [20| , one can conclude that 
this set is sufficient to produce the correct angular dis- 
tribution for (d, p) reactions within FR-ADWA for this 
nucleus. For the (d, p) reactions with 40 Ca at the Ed = 20 
and 56 MeV, DOM performs as well as CH89+WS and 
no significant difference is observed between the standard 
and the dispersive approach (not shown here). 
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There were two cases for which we did see consider- 
able difference in the angular distributions between DOM 
and CH89, the 132 Sn and 208 Pb. These distributions are 
shown in Figs. |5] and [7] It is clear that this is not due 
to the modification in the shape of the overlap functions, 
but rather the optical potentials. The real part of the 
optical potential in the DOM has a larger radius, which 
shifts the diffraction pattern toward smaller angles. How- 
ever an increase in radius in the overlap function only 
affects the magnitude at forward angles (shown by com- 
paring DOM+WS and DOM). While for 132 Sn DOM ap- 
pears to improve the angular distribution, for 208 Pb the 
angular distributions using DOM optical potentials move 
away from the data. We will address the reactions with 
208 Pb in more detail later. 



C. Spectroscopic factors 

While the angular distributions predicted using DOM 
do not differ considerably from those using CH89, the 
normalization of the cross sections do. We determine an 
experimental spectroscopic factor by taking the ratio of 
da I d£L (exp) over da/dCl (theory) for 8 at the first peak of 
the distribution for all but sub-barrier energies. At such 
energies the normalization is determined at backward an- 
gles. We compare the results obtained in the various ap- 
proaches in Table HU the labels CH89+WS DOM +WS 
and DOM are exactly as described in Sec. IIIIBI Note 
that the normalizations are calculated with an overlap 
function normalized to unity. However, the DOM pre- 
dicts overlap functions that are not normalized to unity, 
as correlations are already taken into account. The spec- 
troscopic factor coming directly from such an analysis is 
given m the last column of TableHQand labeled DOM(th). 

One should keep in mind that for 40 Ca, it is notoriously 
difficult to describe low-energy scattering. This fact led 
to the exclusion of proton elastic scattering data at ener- 



TABLE II. Spectroscopic factors obtained from the FR- 
ADWA analysis. The deuteron kinetic energy Ed (lab. frame) 
is in MeV. Reference to the experimental data set used in the 
extraction is also given 



Nucleus 


E d 


data 


CH89+WS DOM+WS DOM DOM(th) 


41 Ca 


20 


[29] 


0.96 


0.85 


0.86 


0.75 




56 


[30] 


0.88 


0.73 


0.74 






2 


W 


0.94 


0.72 


0.66 




49 Ca 


13 


[32] 


0.82 


0.67 


0.61 


0.80 




19.3 


[33] 


0.77 


0.68 


0.62 






56 


[34J 


1.1 


0.70 


0.62 




133 Sn 


9.46 


m 


1.1 


1.0 


0.72 


0.80 


209 pb 


8 


[35] 


1.7 


1.5 


1.2 


0.76 




20 


[36] 


0.89 


0.61 


0.51 





gies below 18 MeV from the DOM fit, as discussed above. 
Since in the ADWA the optical potentials are evaluated 
at E d /2, the DOM results for 40 Ca(d,p) at E d = 20 
MeV are not well constrained by elastic nucleon scatter- 
ing data. In all other cases, DOM potentials provide a 
good description of scattering observables. 

If one focuses first on the traditional CH89+WS ap- 
proach, one sees that a wide range of spectroscopic fac- 
tors can be obtained depending on the beam energy. For 
example with 48 Ca, the spectroscopic factor ranges from 
S = 0.77 to 1.1. This unwanted energy dependence was 
already seen in the systematic study in [28j . Note that 
spectroscopic factors extracted in [28[ do not coincide 
with those in our Table |H] due to additional approxima- 
tions in the reaction theory in [28| , namely the zero-range 
ADWA corrected within the local energy approximation 
(errors introduced in making those approximations have 
been quantified in [l2j ) ■ The large energy dependence is 
significantly reduced when DOM optical potentials are 
used, the exception being the the reactions on 208 Pb, 
which we will consider later. 

When comparing CH89+WS and DOM+WS, one sees 
that the DOM optical potentials reduce the spectroscopic 
factor by a non-negligible amount. Note that this is due 
to the optical potentials alone, since the overlap functions 
arc the same in CH89+WS and DOM+WS. If one now re- 
places the single-particle wave function by that predicted 
from DOM, two things can happen: virtually no change 
(as in the case of 40 Ca) or an important further reduction 
of the spectroscopic factor (as in all other cases) . It is in- 
deed only for 40 Ca that the wave functions are similar. In 
all other cases, the different radius parameter obtained 
in DOM fits, together with the non-locality correction, 
shift density from the interior and enhance the proba- 
bility in the surface region. A larger overlap function at 
the surface, produces larger cross sections which then im- 
ply smaller spectroscopic factors. Values extracted using 
DOM ingredients are much more in line with those from 
(e, e'p) measurements 37j with the exception of 208 Pb. 

Finally, we also show the spectroscopic factors pre- 
dicted directly by the DOM without reference to transfer 
data (last column). These are larger than those extracted 
from experiment for 48 Ca and 132 Sn. This difference is 
mostly associated with the different choice in Ref. [l6| 
for the domain where the imaginary part of the optical 
potentials vanish as compared to the choice in Re f. I I 51) - 
This leads to a larger particle-hole gap in Ref. [16[ as 
compared to Ref. [15| for the same nuclei. Since spectro- 
scopic factors are particularly sensitive to the particle- 
hole ga p, an increase of the order of 0.1 is obtained 
in Ref. [16[ somewhat similar to the difference between 
columns DOM and DOM(th) in Table [ED for these nu- 
clei. If one adopts the perspective that for nuclei where 
optical potentials relevant for transfer reactions are well- 
constrained and the FR-ADWA provides an accurate de- 
scription, the present results suggest that future DOM 
fits should either reduce the particle-hole gap or make the 
coupling at low energy stronger. Nevertheless, the con- 
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TABLE III. The square of the many-body asymptotic nor- 
malization coefficient C 2 for various models (fm _1 ). 



Nucleus E d [MeV] CH89+WS DOM+WS DOM DOM(th) 



41 Ca 


20 


5.0 


4.4 


4.4 


2.8 




56 


4.6 


3.8 


3.8 






2 


31.7 


24.4 


24.4 


29.6 










49 Ca 


13 


27.9 


22.7 


22.6 






19.3 


26.0 


23.1 


23.0 






56 


35.8 


23.5 


23.2 




133 Sn 


9.46 


0.78 


0.71 


0.49 


0.56 


209p b 


8 


4.5 


4.1 


4.2 


2.5 




20 


2.4 


1.7 


1.7 





sistency of the extracted spectroscopic factors employing 
DOM ingredients is very encouraging (see also the dis- 
cussion for 208 Pb below). 

Most of the reactions studied in this work are not 
sensitive to details of the overlap function in the in- 
terior, but rather to the ANC [5|. For completeness, 
in Table IIIII we present the square of the many-body 
ANCs obtained directly from the spectroscopic factors 
presented in Table HH multiplied by the square of the 
single-particle ANC (C 2 = Sb 2 ). The ANC for 41 Ca 
was determined to be Cj7/ 2 = 8.36 ± 0.42 fin -1 , using 
a DWBA analysis of sub-Coulomb (d, p) data The 
analysis in Ref. Q shows that with that large ANC, 
there is no consistency between the expected spectro- 
scopic factor for 41 Ca and the 40 Ca((i,p) 41 Ca data above 
the Coulomb barrier. In Ref. Q it is also shown that 
to obtain consistency between sub-Coulomb and above 
Coulomb energies one needs an ANC about half that 
value. This work solves the discrepancy seen for 41 Ca in 
Ref. @. The ANC for 49 Ca was determined in Ref. [H| 
from a DWBA analysis of sub-Coulomb (d,p) data to 
be C 2 Q/ „ = 32.1 ± 3.2 fm _1 . This value is consistent 
with the values we obtain using CH89+WS. Using the 
DOM instead, the ANC in Ref. [H] should be signifi- 
cantly reduced. An analysis of the 132 Sn(<i,p) 133 Sn reac- 
tion around the Coulomb barrier [39| using FR-ADWA 
and the CH89 optical potentials generates an ANC for 
the ground state of C 2 f7/2 = 0.82 ± 0.07 fin -1 . This 
value would significantly be reduced if DOM were used 
instead. Finally, in Ref. [J] an ANC for 209 Pb is extracted 
from heavy ion reactions C g g/2 = 2.15 ± 0.16 fm^ 1 us- 
ing DWBA. This value is close to the DOM theoretical 
prediction. 

The transfer results for 208 Pb are problematic for two 
reasons: i) there is a very large discrepancy between 
the spectroscopic factors obtained at Ed = 20 MeV and 
Ed = 8 MeV, and ii) for the sub-Coulomb barrier energy, 
this value is much larger than unity. No issues arose 
in the DOM fitting for this nucleus and therefore the 
DOM fits can be considered reliable for the elastic and 
total cross-section data. Even though the transfer an- 



gular distributions are well described within FR-ADWA, 
we cannot exclude the possibility of target excitation. In 
order to test whether the cause for inconsistency is the 
reaction mechanism, we performed exploratory calcula- 
tions including the low-lying 3~ and 2 + states in 208 Pb 
and the strong octupole and quadrupole couplings be- 
tween these states and the ground state within a coupled- 
channel Born approximation (CCBA), using a global op- 
tical potential for the deuteron [|o[. One should keep in 
mind that, within a CCBA approach, the extraction of 
the spectroscopic factor for a single orbital is obscured 
by other components. We found the effects to be strong 
for 208 Pb(d,p) at 20 MeV but rather weak at the lower 
energy We also performed coupled reaction channel cal- 
culations, iterating the transfer coupling, but again found 
the effects to be weak at the lower energy while signifi- 
cant at higher energy. These results suggest that at least 
for the higher energy, the reaction mechanism goes be- 
yond the present implementation of ADWA. The unre- 
alistic SF obtained at the lower energy could indicate 
the failure of ADWA. At present, for reactions involving 
such large Coulomb fields, we cannot verify the validity 
of ADWA through the comparison with exact Faddccv 
calculations. The present techniques used in solving the 
Faddeev equations are limited to targets with Z w 20. 

If the reaction mechanism for 208 Pb(c?,p) at 20 MeV 
necessitates the inclusion of target excitation, could this 
mechanism also play a role in our other test cases, even if 
the angular distributions are well described in our present 
formulation? In Ref. Q a study of target excitation in 
40 Ca(<i,p) 41 Ca at low energy shows strong effects. How- 
ever we expect these to be small at 56 MeV. The 48 Ca 
nucleus has a very weak transition to its first excited 
state, and therefore no significant effect of target excita- 
tion is expected. Target excitation was tested for the re- 
action on 132 Sn and found to be negligible. Nevertheless, 
the present results for 208 Pb do call for an extension of 
the FR-ADWA to include target excitation and deuteron 
breakup in a consistent framework. 



IV. CONCLUSIONS 

We have tested the performance of the dispersive 
optical-model potentials and corresponding overlap func- 
tions in (d,p) reactions. We performed finite-range adia- 
batic calculations for a range of closed-shell nuclei cover- 
ing a wide range of beam energies. Within this descrip- 
tion, the DOM provides all necessary ingredients apart 
from the deuteron V np interaction, i.e. nucleoli optical- 
potentials and the mean-field binding the neutron in the 
final state. We compare the results obtained with the 
DOM to those obtained using a global parameterization 
of the nucleon optical-potential and a standard geometry 
for the neutron single-particle overlap. We find that the 
DOM performs as well as the CH89 parameterization in 
the description of the angular distributions. While spec- 
troscopic factors extracted within the standard approach 
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can be strongly dependent on the energy at which the 
(d,p) data were obtained, this dependence is strongly 
suppressed when DOM potentials arc employed. The ex- 
ception is 208 Pb, a case for which coupled-channel calcu- 
lations demonstrate the importance of target excitation. 
Overall, we find that the extracted spectroscopic factors 
using DOM are significantly reduced when compared to 
the standard approach, bringing the values more in line 
with those obtained from (e, e'p) measurements. 

Because DOM potentials can be extrapolated to rare 
isotopes and separate checks of their quality can be made 
by performing elastic nucleon (proton) scattering experi- 
ments in inverse kinematics, the present framework of the 
DOM in which reaction and structure data are incorpo- 
rated on the same footing, provides an excellent platform 



to analyze transfer reactions involving rare isotopes. 
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